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Abstract 

o\ ■ 

We present a calculation of /b in the static limit, obtained by numerical simulation 
of quenched QCD, at /3 = 6.2 on a 18 3 x 64 lattice, using the SW-Clover quark ac- 
^ ■ tion. The decay constant has been extracted by studying heavy(static)-light correlation 

| functions of different smeared operators, on a sample of 220 gauge field configurations. 

■ We have obtained f^ atlc = (290 ± 15 ± 45) MeV, where the first error comes from the 

. uncertainty in the determination of the matrix element and the second comes from the 

uncertainty in the lattice spacing. We also obtain Mb 3 — Ms d = (70 ± 10) MeV and 
f§®/f§ = 1.11(3). A comparison of our results with other calculations of the same 
quantity is made. 

In the last few years, several groups have computed on the lattice the pseudoscalar de- 
cay constant in the static limit, f B tat , using the method originally proposed by E. Eichten]T0| . 



One can combine the information obtained from the calculation of the decay constant in a 
range of masses near the charm mass with f^ at . In this way it is possible to reduce the 
uncertainty in the prediction of which, on current lattices, can only be obtained by 
interpolation in the heavy quark massf^, [7], [11]. To improve the accuracy on fg it is then 
important to reduce the uncertainties in the static case. These uncertainties are still rather 
large, as it can be seen from the large spread of values reported in the literature[|l|-|| . 
They are mostly due to lattice artefacts, to the calibration of the lattice spacing and to the 
renormalization of the lattice axial current. 

To reduce lattice artefacts, two groups [0, |] have found it convenient to use the SW- 



Clover improved action|12[. In a previous study we have calculated f§ at at (3 = 6.0, using 
cubic smeared sources|§. In this letter we extend the study to a smaller lattice spacing, 
a, corresponding to j3 = 6.2 (again using the SW-Clover action). In the future, a further 
reduction of discretization errors will eventually be obtained by comparing results at different 
values of a and extrapolating to the continuum. This demands a very high accuracy in the 



determination of fp at at any fixed value of the lattice spacing. Our work is a step in this 
direction. 

In reffl, we found that, by using a cubic smearing source in the Coulomb gauge, there 
is an optimal smearing size for which it is easiest to isolate the lightest pseudoscalar state, 
corresponding to the B- meson in the static limit. We confirm this result and find that, 
using a double smeared source[[L], § at (3 — 6.2, the optimal size is L s = 11 (L s = 13 is the 
optimal size for single smeared sources). L s = 11 and double smeared sources have been used 
to obtain the results reported in the abstract and belowQ We will also compare different 
smearing sizes (7 < L s < 13). 

The decay constant in the static limit can be derived from the matrix element of the 
static-light quark lattice axial current Z L 0: 



f s P tat sjM P = V2 Z s 2 at {M P a, a s ) Z L a" 3/2 (1) 

where fp tat is the decay constant for a meson of mass Mp. The renormalization constant 
Z^ at (Mpa,a s ) includes the anomalous dimension of the axial current in the effective the- 
ory and the 0(a s ) correction necessary to relate the lattice operator in the effective theory 
to the continuum operator in the full theory. The estimate of Z^ at is subject to a large 
uncertainty ||T3||- |T^| . To evaluate it, we have used the "boosted", tadpole improved per- 



turbation theory proposed in ref. |18], |19~ . The uncertainty in the determination of Z^ at is 
given by the uncertainty on the value of the boosted coupling g v to be used and by higher 
order corrections. We have used two different definitions of g v , i.e. gy = (8i^ c ) 4 6//5 and 
9v = Q/p/ < l/3TrU p >§, H, 0. With the SW-Clover action, unlike the Wilson case, the 



two determinations of g v almost coincide and we get Z^ at = 0.805 — 0.812 [1151 [TBI]. Thus we 
have chosen Z^ at = 0.81. Using a^ 1 = (3.0 ± 0.3) GeV, as suggested by the calibration of 
the scale from M p and from this work and ref.[^8]|, together with Z L = 0.111(6), we get 



f s P tat y/Mp~ = (0.66 ± 0.04 ± 0.10) GeV 3/2 corresponding to: 

fg at = (290 ± 15 ± 45)MeV, (2) 

where the first error comes from the uncertainty in Z L and the second comes from the 
uncertainty in the scale. When the scale is taken from the string tension, a -1 ~ 2.7 GeV, 
one gets a lower value, i.e. f§ at = (245±15) MeV. The calibration of the scale from the string 
tension is more convenient when studying the scaling properties of fp at with a, since the 
error on this quantity is of 0(a 2 ). It is however difficult to express the lattice string tension 
in terms of a well defined experimental quantity, in order to obtain the absolute scale. If we 
use the naive perturbative value of the renormalization constant, ie. Z^ at = 0.88, and use 
a- 1 = (3.0 ± 0.3) GeV, we obtain instead f s B tat = (315 ± 15 ± 45) MeV. 

Most of the methods and techniques used in the present study have been developed and 
explained in the literature on the same subject [|T|-0. Details can be found for example in a 
previous publication by our collaboration ||. In the following we explain only those aspects 
which are new with respect to ref.[Q. The main differences are the following: 

1) smearing in the origin; 

2) calculation of the matrix elements of the local axial current, Z^, by fitting the smeared- 
local correlation function; 



1 At (3 = 6.0 the optimal smearing size was found to be L s = 5 — 7. For a more detailed discussion see 
ref.§. 

2 The precise definition of Z L will be given in eq.(||). 



3) a method to extract Z^ from ratios of correlation functions and without any fitting 
procedure. 

In this study, we have used local-smeared correlation functions with the smeared source 



in the origin. This strongly reduces the statistical noise in the correlators [20]. In our previous 
work at (3 = 6.0 we used instead the smeared-local correlations with the smeared source in 
the sink. To improve the results of ref. ||, we are planning to repeat the calculation at 
(3 = 6.0 with the smeared source in the origin. Points 2) and 3) will be discussed below. 

We have generated 220 configurations in the quenched approximation, at — 6.2 on 
a 18 3 x 64 lattice. The configurations were obtained from several independent runs, each 
with a thermalization of 3000 Montecarlo sweeps (5 hits), starting from a cold configuration. 
Independent configurations were gathered every 800 sweeps. We have then transformed each 
configuration to the Coulomb gauge, with the same accuracy as in ref. [BJ . We have used the 



0(a) improved SW-Clover action|T^]. The operators appearing in the correlation functions 



are "improved" by modifying the light quark propagators as explained in sec. 9 of ref.Q. 
The statistical errors quoted in this letter have been estimated using the jacknife method, 
by decimating 10 configurations at a time. In ref. || we verified that this gives a reliable 
estimate of the errors. Throughout this letter we have computed light quark correlation 
functions using a "thinning" trick. This method consists of the following: one only considers 
one point out of three in each direction when summing the correlation functions over space, 
at fixed time. This allows a saving of a factor of 27 in the memory required for the quark 
propagators. Thinning was a necessity imposed by memory limitations of the APE-tubo[pi[] 
used for the present simulation. We have checked that the effect of thinning on 2-point 
correlation functions is totally negligible | 



We have inverted the light quark propagators for three values of the bare quark masses, 
K = 0.14144,0.14190,0.14244. For comparison, the first K value has been chosen equal to 
one of those used by the UKQCD collaboration^, at the same value of (3. Our lightest 
quark mass is slightly larger than in ref.Q since, for lack of memory, we have to work on 
a smaller space volume ( the UKQCD Collaboration is working at /3 = 6.2 on a 24 3 x 48 
lattice). 

We have measured the pseudoscalar (M 5 ) and vector (My) meson masses and the matrix 
elements \/~Zl =< 0\^ 5 ip\P > and y/Zy =< 0\ijj'yk' l l , \V > by fitting the correlation functions 
of the local pseudoscalar density and local vector current with components k — 1, 2, 3, in 
the time interval 15 < t < 28. By studying several time intervals, we have checked that the 
contamination due to higher mass excitations induces a negligible systematic error on the 
masses and the matrix elements of the local operators, in the range of quark masses explored 
in this study. The results at each K and K c ( the critical value of the hopping parameter) 
are reported in table |I[ From a linear fit in 1/K to Mf we obtain K c = 0.14315(4). Using 
Mk = 498 MeV and the scale from the rho mass we obtain K s = 0.1422(2), where K s is the 
Wilson parameter corresponding to the strange quark. 

The physical size of the lattice, obtained from the vector meson mass extrapolated lin- 
early in Mf to the chiral limit, is a -1 = (3.2 ± 0.3) GeV. The pseudoscalar (/p) and vector 
(1/ fy) decay constants have been extracted using standard methods, see for example ref.||, 
from the local vector and axial vector currents. In table |l| we report the lattice values 
of fp/Z^ en , the ratio fp/Z^ en My and 1/Zy en fy as a function of the light quark mass. 
Zy e £ are the renormalization constants of the vector and axial vector lattice currents [^3 - 
which, for the SW-Clover action, have been computed non-perturbatively at /3 — 6.2, 



K 


0.14144 


0.14190 


0.14244 




M 5 


0.292(3) 


0.247(3) 


0.187(4) 




My 


0.375(6) 


0.338(8) 


0.291(18) 


0.238(21) 


Z 5 x 10 3 


8.55(44) 


7.43(46) 


6.65(60) 




Z v x 10 3 


1.80(18) 


1.28(18) 


0.76(27) 




fp/Zf n 


0.063(2) 


0.059(2) 


0.052(3) 


0.045(4) 


fp/Zf n M v 


0.169(5) 


0.173(7) 


0.180(15) 


0.186(20) 


1/Z«* n fv 


0.302(7) 


0.313(10) 


0.326(22) 


0.341(27) 



Table 1: Pseudoscalar and vector meson masses and matrix elements at different values of 
K . The corresponding values extrapolated linearly to K c are also given. 



Z*} en = 0.831(2) and Zf n = 1.04(1) [[?§]. We have extrapolated f P /Zf n , f P /Zf n M v and 
1/Zy en fy linearly in Mf to the chiral limit, obtaining the numbers in table [l| Using the 
non-perturbative determinations of Zy en and Z^ en given above, and M p = 770 MeV we get 

U = (149 ± 16) MeV , — = 0.28 ± 0.02 (3) 

Jp 

We also obtain f K /fn - 1 = 0.11 ± 0.04 and 1/% = 0.27 ± 0.02. We give for comparison 
the experimental values: / w = 132 MeV, \ jf p = 0.28, f K /f v - 1 = 0.22 and l// = 0.23. 
The slope of the pseudoscalar decay constant as a function of the quark mass is lower than 
what we found in previous calculations at (3 — 6.0 and for this reason we get a rather low 
value for fx/ ' f n — l. We believe however that this effect will disappear with better statistics. 
We can also use f n to fix the scale. By using the value of Z^ en quoted above we obtain 
a -1 = (2.8 ± 0.3) GeV to be compared to the value obtained from M p of a' 1 = (3.2 ± 0.3) 
GeVQ. At /3 — 6.0, we found that the scale determined from f n and M p were in very good 
agreement. In the present case there is a difference of about 400 MeV that we cannot resolve, 
given the large statistical errors. We also obtain M* K = (893 ±7) MeV and M = (1020 ±13) 
MeV which compare very well with the experimental values of 892 MeV and 1019 MeV 
respectively. 

The static B meson decay constant, f^ at has been obtained from the correlation functions 
of three different currents: 

A L p (x) = Q(x) lpl5 q(x) (4) 
A S p (x) = J2Q( x ih^q(x) (5) 

i 

A°(x)=Y,Q( x i)wi( x i) (6) 

where the sum over % (J) is extended over a spatial cube of size L s , centered at x and the 
heavy quark Q is taken from the static limit. We have computed the correlation functions 
of the above currents, C LL } C LS , C LD , C ss and C DD } for L s = 7,9, 11, 13. C LD is defined 
as: 

3 From the axial current-axial current two point function one can also get a lattice estimate of / w . In this 
case, using the same Z^ en and the experimental value of /„. we also get a -1 = (2.8 ± 0.3) GeV. 



c LD (t) = £< 0|^(f,tX(0,0)t|0 > (7) 

X 

and similarly for the others. At large time distances the correlation functions behave as the 
propagator of a single state: 

C LD (t) = Z L Z D e- AEt (8) 

The constants Z L , Z s and Z D are related to the matrix elements of the different axial currents 
and AE is the binding energy of the lightest pseudoscalar state. 

A standard method to extract Z L is the following: 

1) At large time distances the correlation functions are dominated by the lightest prop- 
agating state. Let us call ti the time at which we start observing a plateau both for the 
effective binding energy: 

AEgi = log[C DD (t)/C DD (t+l)] (9) 

and for the ratio: 

ti clearly depends on the smearing size and it is determined by observing that AE DD and R 
do not vary appreciably from time to time. 

2) For t > ti, i.e. in the time interval where (|5p and (ffD) have a plateau, we fit the 
two-point smeared-smeared correlation C DD to the expression: 

C DD (t) = (£ D )V AW (11) 

From this fit we obtain the matrix element of the corresponding axial current, Z D . 

3) We then compute R defined as: 

R _ gg* m/* 2 m (12) 

where <jR(t) is the jacknife error on the ratio R(t) at the time t. We used tf = 15 for which 
the errors in the static-light correlation functions are reasonable. 

4) Z L is then obtained from the product: 

Z L = Rx Z D (13) 

where Z D is found from the fit of the correlation function to the expression in eq. (|TT]) . We 
will call this method DD-mass fit. 

We now give two other methods that we have used in this letter to determine Z L . 

i) LD-mass fit: in point 2) we replace the C DD fit with a fit of C LD . In this way we 
obtain the binding energy and P = Z L Z D , i.e. the products of the matrix elements of the 
local and smeared currents. We then compute Z L from the product of P with R as defined 
in 3), Z L = y/R x P. The advantage of this method is that C LD has the smaller statistical 
errors. A possible disadvantage is that the effective binding energy for C LD approaches its 
asymptotic value at large times from below, contrary to the case in which one uses C DD . 
However the quality of the plateau that we observe in the effective mass AE LD , from C LD , 
is such that we do not doubt that we have isolated the lightest state. This is shown in 



Figure 1: AE^ and AE^f) + 0.2 from C DD and C LD as a function of the time t for L s — 11 
and K = 0.14144. We also report the published data of the UKQCD collaboration^ . The 
points of ref. ^ have been moved by —0.2 in order to distinguish them from our results. The 
lines joining the points are there to help the reader distinguish the different cases. 



fig.[l] where AE LD and AE DD are shown together (see also table 0). For comparison on the 
quality of our data we also report the results of the UKQCD collaboration, at the same value 
of (3 and K } obtained with a statistics of 20 configurations on a 24 3 x 48 lattice. We will 



compare our static results with those of ref. [0] and our light-light results with ref. [£§] at the 
end of this letter. 

ii) Ratio method: we compute the following ratio of correlation functions: 

iWMa)- C ™(t 1 + t 2 ) (M) 

At large t\ and t 2 (ti,2 > U), it is clear from eqs.(§,[Ll]) that R z l — > (Z L ) 2 . In practice, we 
average R z i,(ti, t 2 ) from ti + t 2 = t m = 2ti to a certain maximum value t± + t 2 = tu-, for which 
the errors are reasonably small, with all ti j2 which satisfy t\ i2 > ti. The advantage of this 
method is that Z L is obtained directly and no fitting is necessary. The possible disadvantage 
is that we have to compute C DD at large time distance, i.e. > 2t,. We note here that this 
method can be applied not only to measurements of fp at , but to a large variety of matrix 
elements. 

Several consistency checks are possible. One should find a plateau corresponding to the 
same value of the effective mass for the C LD and C DD . Since we extract Z L by taking ratios, 
we also need to overlap the time ranges in which the plateau in the effective mass is observed. 
We show that this is possible in figs.[l],|2] where we give AE e Jl and AEffi function of 
the time for L s = 7 — 13. 



Figure 2: AEpD and AEfy + 0.2 from C DD and C LD as a function of the time t for 
L s = 7 - 13 and K = 0.14144. We also give best estimate, AE eff = 0.561(3), the band of 
which appears as a thick line in the figures. Notice that L s — 11 works well both for C LD 
and C DD at moderate values oft. 



A more quantitative check of the reliability of the results is given by the fact that different 
correlation functions give the same effective binding energy. The results and errors for AE e Jl 
and AE^ with different smearing sizes and time intervals, which we consider in the good 
range of L s and t, are reported in table 0. The binding energies extrapolated linearly in M| 
to the chiral limit are also given in the same table. From the dependence of the binding 
energy on the quark mass one can derive the mass difference M Bs — Me d |29|]. Our best 
estimate for this quantity is M Bs a - M Bd a = 0.020(4), ie. M Ba - M Bd = 70(10) MeV. We 
can see from fig.0 and table |2] that the smearing with the best plateau both in the DD and 
LD case is L s = 11. Consequently all the results given for fp attc have been obtained with 
this smearing size and ti = 5 and 7. 

A further quantitative check is that the three different methods give the same value of 
Z L . A comparison of the different determinations can be done from table [| Z L , extrapolated 
linearly in M| to the chiral limit, is also given in the same table. From the dependence of 
Z L on the quark mass we have obtained the value of f B J f Bd = 1.11(3) (see 0). Finally, 
we have checked that the S— smearing (see eq.(|)), when we use L s = 13, gives consistent 
results to the D— smearing (eq.(|[)), thus giving confidence in our methods. 

From the results for the effective mass and Z L at different values of the Wilson parameter 
and in the chiral limit, other groups working on the same problem and using the same quark 
action, but different smearing techniques to isolate the lightest state, can check our results. 
Our best estimate of Z L is 0.111(6) which we obtain from a weighted average of the results 
in table |3], where the error includes the spread of values in the table added in quadrature 
with the statistical error. 





AEdd 


K 


L a =Ut = 5-15 


L s = llt = 7-15 


L s = 13 t = 5-15 


L s = 13 t = 7-15 


0.14144 
0.14190 
0.14244 


0.560(5) 
0.548(6) 
0.532(6) 


0.556(7) 
0.545(7) 
0.529(8) 


0.563(6) 
0.551(7) 
0.535(8) 


0.560(8) 
0.548(8) 
0.531(10) 


K c 


0.514(6) 


0.511(9) 


0.515(7) 


0.513(10) 




AE LD 


K 


L a =Ut = 5-15 


L 8 = \lt= 7-15 


L s = 13 t = 5-15 


L s = 13 t = 7-15 


0.14144 
0.14190 
0.14244 


0.562(2) 
0.552(3) 
0.541(3) 


0.561(3) 
0.551(3) 
0.539(3) 


0.550(2) 
0.540(3) 
0.528(3) 


0.553(4) 
0.543(4) 
0.531(5) 


K c 


0.526(4) 


0.524(5) 







Table 2: AEdd and AEm for L s = 11 and 13. The fits done in different time intervals 
have been reported in order to show the stability of the results. The agreement is very good 
with the only exception of AEld when L s = 13. The reason can be understood by looking at 
fig$ where it is shown that there is a drift in the effective mass as a function of the time. 
Analogous problems are encountered with C LD for L s = 7 and 9. On the contrary this table 
shows that with L s — 11 one gets results consistent, within very tiny errors, from both C LD 
and C DD and different time intervals. 



Once Z L is established for a given value of the 13, the determination of f^ at still depends 
upon several factors. Apart from effects of 0(a), which we expect to be smaller with the 
SW-Clover action, there are the uncertainties coming from the calibration of the lattice 
spacing, a, and from the renormalization constant of the axial current in the static theory 
(see eq. 

There is only one group which has worked at this value of (3, using the SW-Clover action, 
UKQCD|7], |2"8f . We now compare our results with theirs for the light meson sector and the 
static case. 

In figs.| we plot Mf and My clS Sl function of l/K. Our results are in good agreement 
for Mf (we find essentially the same K c ) but are systematically lower for My. In fact, at 
the only value of K which is common to the two calculations, our result coincides with the 
value of My reported as My w in table 1 of ref . [J28| . The results are in better agreement if 



we include in the comparison the systematic error quoted in ref. ||28|| , which was estimated 
by looking at the different determinations of My, obtained by varying the time interval of 
the fits. A possible origin of the difference is that UKQCD fits the correlation functions 
from ti = 13 to 23. The difference in My is reflected in the estimates of a -1 , the inverse 
lattice spacing, which in their case is a -1 = (2.7 ±0.1) GeV (no systematic error is reported) 
to be compared with our value of a~ 1 = (3.2 ± 0.3) GeV. On the contrary, our result for 
fp is systematically above the result of ref.[^], see also table [I] and figs. [| However the 
results are in this case compatible within the statistical fluctuations. The combined effect 
of a smaller My and a larger fp is that the difference between the two groups becomes a 
discrepancy for fp/Myf\. From the above discussion we have seen that most of this is due 

4 Note that we find that fp/My/Z^ en is slightly increasing in the chiral limit, contrary to previous findings 
at P = 6. OH. This is probably due to a statistical effect. 



Z from the DD-mass fit 



K 


L a = \lt = 5-15 


L a = llt = 7-15 


0.14144 
0.14190 
0.14244 


0.134(3) 
0.127(3) 
0.119(3) 


0.132(4) 
0.125(4) 
0.116(4) 


K c 


0.108(4) 


0.104(5) 




Z L from the LD-mass fit 


K 


L s = Ut = 5-15 


L s = 11 t = 7-15 


0.14144 
0.14190 
0.14244 


0.136(2) 
0.130(2) 
0.123(2) 


0.135(2) 
0.129(2) 
0.121(2) 


K c 


0.114(2) 


0.111(3) 




Z L from the Ratio method 


K 


L s = 11 t = 10-15 


L s = 11 t = 14-15 


0.14144 
0.14190 
0.14244 


0.135(2) 
0.129(2) 
0.121(2) 


0.134(5) 
0.128(5) 
0.119(6) 


K c 


0.111(2) 


0.109(7) 



Table 3: Z L obtained from the three methods (see text) using L s = 11. The fits are done in 
two time intervals (t — 5 — 15 and t = 7 — 15) in order to show the stability of the results. In 
the ratio method we have chosen t m = 10 and tu = 15 for ti — 5 or t m — 14 and tu = 15 
for ti = 7, see text. The agreement between all three methods is quite satisfactory. 



Figure 3: Mf, My, f P /Zf n and f P /(Zf n M v ) as a function ofl/K. We plot both our and 
the UKQCD^ZBJ results, the latter obtained with a statistics of 60 configurations. Only the 
statistical errors are shown. The lines are linear fits in 1/K to our data. 



to Mi 



v- 



For the static case, at K=0. 14144, our results for the binding energy are definitely 
incompatible with those of UKQCD]?]]. They quote AE e ^ = 0.59 ± 0.1 to be compared 
with our results, given in table |2|. At this point a comparison of Z L becomes difficult. In 
any case, the difference for Z L is not too large: they find Z L = 0.142lg at K=0. 14144, to be 
compared with a typical value (error) for our results of Z L = 0.136(2), and Z L = 0.1241? at 
K c , to be compared to our estimate of Z L = 0.111(6). 

A prerequisite to a precise determination of f^ at demands agreement on the value of 
AE e ff and Z L obtained with a given action at a certain value of (3. Our group and the 
UKQCD collaboration have decide to work with the SW-Clover action in order to reduce 
lattice artefacts. In this paper we have presented a very accurate determination of Z L , 
derived with three different methods, which give compatible results, within tiny statistical 
errors. Our results for the binding energy show some difference with the results of ref.f/J 
which, by taking the quoted statistical errors, appear to be of systematic origin. The same 



can be said for the mass of the vector meson p8| . Other differences, for example fp, can 
instead be the effect of statistical fluctuations. 
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